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Abstract 

We give a general formula for the cross section for diffractive vector meson electro- 
production, 7*p — ► Vp. We first calculate diffractive qq production, and then use parton- 
5_i ■ hadron duality by projecting out the J p = V~ state in the appropriate mass interval. 



We compare the Q 2 dependence of the cross section for the diffractive production of p 
and J/ip mesons with recent HERA data. We include the characteristic Q 2 dependence 
associated with the use of the skewed gluon distribution. We give predictions for ctl/ot 
for both p and J ftp production. 



The diffractive leptoproduction of vector mesons at high energy is an interesting and im- 
portant process. Indeed diffractive 7*p — > Vp data, with V = p,u,<p, J /ip and T, are becoming 
available with increasing precision from the experiments at HERA They offer the op- 

portunity to study the vacuum-exchange singularity as a function of the virtuality Q 2 of the 
incoming photon and of the mass M of the produced vector meson. Moreover observation of 
the vector meson decays allows both and o>r to be measured, and even s-channel helicity 
conservation to be checked ]|, || ]8| , (| . 

Let us first review the description of diffractive electroproduction of p mesons; a pro- 
cess which has attracted a lot of theoretical interest JTU|-[I1|. At first, phenomenological 
parametrizations based on the vector-meson-dominance model and Regge exchanges were used. 
Then a non-perturbative two-gluon exchange model of the Pomeron was introduced flTUf . For 
large Q 2 however, we would expect that a pure perturbative QCD description is applicable. 
Such a description for the production of longitudinally polarised p mesons was given by Brod- 
sky et al. |IT| , using the leading twist wave function for the p meson. The process is sketched 
in Fig. 1. However for the production of transversely polarized p mesons, the perturbative 
QCD approach encounters an infrared divergence in the integration over the quark transverse 
momentum. This problem can be overcome by using parton-hadron duality [Tj|. The wave 
function of the p meson then never enters explicitly. The only property that is used is that 
the p meson corresponds to the J p = l~ projection of 'open' qq production (with q = u,d). 
The projection has the effect of curing the infrared divergence. The resulting cross section is 




Figure 1: Schematic diagram for diffractive vector meson production at HERA, 7*p — > Vp. 
The longitudinal fractions x and x' of the ingoing and outgoing proton momentum carried by 
the gluons are given by Eq. fllCf); the gluons have momenta ±£t transverse to the proton, z 
and 1 — z are the longitudinal fractions of the photon momentum carried by the q and q, and 
±kr are their momenta transverse to the photon. There are four possible couplings of the two 
gluons to the q and q, represented by the upper circle. 



then integrated over an appropriate interval AM of the invariant mass of the qq pair which 
covers the p resonance peak. As there are almost no other possibilities^ for hadronization of 
the qq pairs at M q q ~ M p , the procedure is expected to give a reasonable estimate of the cross 
section for p electroproduction. Indeed this perturbative framework |13[ was found to describe 

x We allow for ui production by taking the ratio w : p to be 1:9. 
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the Q 2 dependence of p electroproduction for Q 2 > 5 GeV 2 for both longitudinally and trans- 
versely polarised p mesons, including the observed Q 2 dependence of the (Jl/gt ratio. The Q 2 
behaviour of the amplitude is governed by the structure of the quark propagators and by the 
effective anomalous dimension 7 of the gluon, defined by xg(x, K 2 ) ~ (K 2 ) 1 . In particular the 
naive expectation that = o^M 2 jQ 2 is modified to0 




(1) 



which, on account of the decrease of 7 with increasing Q 2 , is in good agreement with the 
observed o L fo T behaviour with Q 2 . 

The cross section for the diffractive electroproduction of vector mesons is proportional to 
the square of the off-diagonal or skewed gluon distribution. That is x 7^ x' in Fig. 1, whereas 
for the conventional (diagonal) gluon we have x = x' . In fact the gluon distribution becomes 
more skewed as Q 2 increases, and is more skewed for vector mesons of larger mass M. Skewed 
distributions were not included in our predictions of the Q 2 behaviour of p electroproduction in 



Ref. [|13|]. When compared with subsequent precise HERA data ||, these predictions were found 
to fall off a bit too rapidly with increasing Q 2 . We will see that the effect of using the skewed 
distribution, rather than the usual approximation of using the conventional (diagonal) gluon, 
will enhance the cross sections at the larger values of Q 2 at which data exists. As was previously 



discussed [p~3| , there are uncertainties in the normalization of the predictions of the diffractive 
cross sections, but much less in the predictions of the energy or Q 2 dependence. Nevertheless, 
in order to use the Q 2 dependence of the data to reveal the effects of the skewed distribution, we 
must include the Q 2 dependence of other effects in our calculation. The (imaginary part of the) 
amplitude is calculated at t — and the cross section obtained by integrating da/dt ~ exp(— bt) 
over t. We must therefore allow for the decrease of b with increasing Q 2 . Second we must 
study the ambiguity in our estimates of the next-to-leading order (NLO) correction. In the 
perturbative region, we find that the Q 2 variation of p electroproduction from these two sources 
is smaller than that due to the use of the skewed gluon distribution. Also we must, of course, 
include the contribution from the real part of the amplitude. When we compare the full QCD 
prediction with the Q 2 behaviour of diffractive p and J/i[) production recently measured at 
HERA we find that the data are compatible with the characteristic enhancement arising from 
the skewed gluon. 

We use perturbative QCD to derive the general formula for the cross section for diffractive 
vector meson production by first recalling the formula for diffractive production of a qq system 
of mass M. For production from a transversely (longitudinally) polarised photon 



d 2 a T ^ 



dM 2 dt 



27r 2 e 2 q a 



t=o 



3(Q 2 + M 2 ) 2 



dz 



B T(L) 2 



(2) 



2 Eq. (Q) is an approximate result obtained assuming that, for each Q 2 , 7 is constant throughout the in- 
tegration over the quark loop. The full calculation can be found in Ref. jl3|]. See also the results discussed 
below. 
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where i = +,— and i' = +, — denote the helicity of the quark and antiquark. The helicity 
amplitudes are 



ImBi 



T 



2jz(l-z) 



B 



ImB 1 



0. 



(3) 



(1 — z)kxlT 



'z(l-z) ^(l-z) 
for a photon of helicity +1, whereas for a longitudinal photon we have 

B^ + = B L _ = 0, 



(4) 



ImBi 



-lmB L 



'z(l - z)Q 2 



It, 



The variable z is the fraction of the photon's momentum carried by the quark, is the 
transverse momentum of the quark relative to the photon, e q is the charge (in units of e) and 
m the mass of the quark; a = 1/137. The mass M of the qq system satisfies 



M 



m 2 + k% 



z(l-z) 

The integrals over the transverse momentum ±£j< of the exchanged gluons are [13, |H| 



(5) 



h(K 2 ) 
h{K 2 ) 



fa s (£ 2 T )f(x,x',f T ) (l - 



K 2 



a s (£ 2 T )f(x,x'J 2 T ) 



K 2 2k 2 T 



+ 



where 



K 



z(l - z)Q 2 + k 2 T + m 2 



(6) 
(7) 

(8) 



K 2 



(K 2 + £ 2 T ) 2 



(9) 



The function f(x, x', £%) is the skewed unintegrated gluon distribution describing the lower part 
of Fig. 1. The momentum fractions carried by the exchanged gluons satisfy ||T5 



x 



Q 2 + M 2 ^ 
W 2 + Q 2 t 

where W is the j*p centre-of-mass energy. 



> x' 



W 2 + Q 2 



(10) 
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In the strict leading log(l/x) approximation, it is enough to use the diagonal unintegrated 
distribution, as at each splitting we keep just the leading log(l/x) terms and neglect the cor- 
rections due to x' <C x. In this limit 



f(x, x , £ 2 



f(x,£ 2 ) 



d(xg(x,/j, 2 )) 



<91n /i 2 



(11) 



and there is no difference between the diagonal and skewed distributions fll6| . This is the 
approximation which is conventionally used. 

Here we wish to take into account the skewed effect, but first we must extend the definition 
of the unintegrated gluon, (|TT|), beyond the leading log(l/x) approximation. Indeed it is easy 
to see that ( |TT] ) can only be true for sufficiently small x. If x increases then / calculated from 
( [TTD would soon become negative due to the (negative) virtual contribution in the DGLAP 



evolution. It was shown in Ref. [|17[ that the virtual corrections may be resummed via the 
Sudakov form factor and that the number of gluons with transverse momentum It is 



f(x,f T ) 

where in the double log approximation (DLA) 



d[xg(x,q%)T(q$,n 2 )] 



d In go 



a 2 =E 2 



(12) 



exp 



-C A a s (fi 2 ) 2 
4tt n ql 



(13) 



with scale fi 2 ~ (Q 2 + M 2 )/4. T is a survival probability. It is the probability that the parent 
gluon does not emit gluons in the interval q^ < q T < fi 2 . From the formal point of view, the 
T factor may be regarded as a next-to-leading order correction since the main contributions to 
the integrals and (|7|) come from the region^ £ T < fi 2 . In general we find that the inclusion 
of T has a small effect, essentially only ensuring the positivity of / for T production which 
samples values of X clS large as x ~ 0.05. 

The main effect of using the skewed (or off-diagonal) gluon distribution comes, within leading 
InQ 2 kinematics, from the region where x' -C x, see (|10|). In this region the skewed gluon 
distribution H g (x,x') (integrated over £?) is larger than the conventional diagonal distribution 
H g (x,x) = xg(x). For small x, which is appropriate for vector meson production at HERA, 
the enhancement is generated entirely by off-diagonal evolution. Moreover the ratio 



Rn 



H g (x, x' x) 

Hg(X,X) 



(14) 



3 If It > H then we set T = 1 in (^2|), consistent with the DLA. It may occasionally happen (at the edge of 
phase space) that the inclusion of the T factor in the DLA is not enough to ensure the positivity of f(x,lj), 
whereas the exact form of the T factor would guarantee that / > 0. Therefore we set / — if it should happen 
that is negative. 
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can be determined unambiguously in terms of the known diagonal distribution [18]. It was 
shown that the enhancement R g depends on the effective power X(Q 2 ) of the small x behaviour 
of the gluon xg ~ x~ x . The result is |H| 

2 2A+3 r (A + f ) 

Note that the off-diagonal enhancement enters at leading order (in InQ 2 ) and increases with 
Q 2 (since A increases with Q 2 ). Here we allow for the off-diagonal effect by multiplying the 
amplitudes (ff) and (f|), calculated^ with diagonal gluons in @ and (0), by the factor R g . We 
determine the effective power for each component amplitude separately, that is 

d\ogB iV 

A " 01og(l/s)- (16) 

The full NLO corrections for the diffractive process are not known yet, and so we approxi- 
mate them by a /C factor |14], [13j . Following [TJj , the main {ir 2 enhanced) part of the K, factor 
is of the order of (7r 2 CVas/7r), where Cf = 4/3. It comes from the iir terms in the double 
logarithmic Sudakov form factor exp [-C F (a 5 /47r)ln 2 (-M 2 )] where ln(-M 2 ) = InM 2 + in. 
Thus we multiply the amplitudes by the factor [|TJJ] 

K = exp(7rC F a s /2). (17) 

But we still have the ambiguity of the choice of the scale of a$- We show results for two choices 
of scale: p 2 = K 2 and our default value fi 2 = 2K 2 , where K 2 is given by (|8|). 

To include the contribution from the real part of the amplitude we use the signature factor 

£(+) = i + tan(7rA/2) (18) 

for positive signature exchange. This is a simple way of implementing the dispersion relation 
result. It gives 

ReB^ = tan(vrA/2) JmB w , (19) 



where A is given by (|16f) . The inclusion of the real part enhances the cross section of p production 
by 14 to 19% in the range where we compare to data, J/tp production by 18 to 25%, and T by 
about 30%, where the bigger effect always occurs at higher Q 2 . 

So far we have calculated the cross section d 2 a/dM 2 dt for diffractive qq production at t — 0. 



To determine da/dM we integrate the form exp(— bt) over £, with [JT9 

2 4 2 / W 2 M 2 \ 

KQ } = ((t) + 0.71 GeV 2 ) + Q 2 + M* + (t) + 2a ' ln { (Q 2 + M^ ) ' (20) 



4 In the infrared region l\ < £q we use the "linear" approximation = (i"^/ ^)as{(-f i )g{x 1 £q) as 

described in This linear approximation is reasonable since (i) it corresponds to a constant gluon-proton 
cross section at small scales It < £o and (ii) it matches well to the scale dependence of the phenomenological 
gluon distribution at low It- We have checked that our results are stable to reasonable variations of £\ about 
our default value of 1.5 GeV 2 . 
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where (t) is the average value of t. (Here we set (t) = 0.) This form, with a'p = 0.15 GeV -2 , 
successfully reproduces the t behaviour of diffractive p meson leptoproduction data as a function 
of Q 2 and M 2 ~ M v . It is motivated by the additive quark model, together with a form factor 
given by Fy{t) = M 2 /{M 2 — t), see also PD[. We also use the phenomenological expression 



(|20f) for diffractive J/ip production even though the measured slopes appear, at present, to be 
about 2 GeV~ 2 or 30% less. Using the observed values of b would lead to an overall increase 
in the J/ip cross section of about 30%, well within the present uncertainties in the theoretical 
normalization. 

To determine the cross section for 7*p — > Vp from that for diffractive qq production, we 
project out the J p = 1~ state in the qq rest frame. However the helicity amplitudes B^i are 
defined in the target proton rest frame, and helicity is not conserved by Lorentz transformations 
for the heavy quark states. So to obtain the helicity amplitudes Ajji in the qq rest frame for 
V = J/ip and T, we must perform a Lorentz boost and use 

A«' = Ci 3 C i' i ' Bii', (21) 



where the known coefficients Cy are given in Ref. pi| . Finally we integrate the cross section 
da/dM 2 for J = 1~ qq production over an appropriate interval AM 2 covering the vector meson 
resonance. Clearly this, together with the K, factor of (|TTD, introduces an overall normalization 
uncertainty. However here we are interested in the Q 2 dependence of cr(7*p — > Vp) and the 
properties of the ratio <Jl/&t, rather than the normalization. 

The QCD predictions for p, J/ip and T production^ are compared with HERA data in 
Figs. 2-4. Fig. 2 shows the predictions obtained from using three different recent gluon dis- 



tributions: MRST [g, CTEQ and KMS All fit the F 2 structure function data well. 
The first two are obtained from conventional NLO DGLAP analyses, while the KMS analysis 
is in terms of an unintegrated gluon distribution which satisfies a unified BFKL/DGLAP equa- 
tion with subleading ln(l/x) contributions. In each case the scale p? in as of the K, factor is 
taken to be 2K 2 . The lower continuous curves in Fig. 2 correspond to the 'default' prediction 
obtained using MRST99 partons. In both plots the upper continuous curve corresponds to 
taking p 2 = K 2 . The dashed curve is the default prediction using the diagonal gluon, and so 
comparison with the lower continuous curve shows the enhancement due to off-diagonal effects. 
At the larger values of Q 2 the enhancement is about 55% for p production and 70% for J/ip 
production. 

Our approach is infrared finite. However there are non-negligible contributions from the 
region of low gluon transverse momenta £t < A)- Fortunately the predictions based on conven- 
tional DGLAP partons (MRST, CTEQ) are rather insensitive to the choice of the value of £q. 
Nevertheless to check the infrared sensitivity of the predictions we also use a gluon obtained 

5 The J/ip production amplitudes were calculated for a charm quark mass m c = 1.4 GeV, whereas for T 
production we take the b quark mass nib = 4.6 GeV. 
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from a unified BFKL/DGLAP analysis of the deep inelastic data [24]. Q We would expect some 
difference since the latter (Reggeised) gluon embodies a higher twist component originating 
from the BFKL evolution, which may be important at low scales. Indeed we see from the 
dotted curves in Fig. 2 that the cross sections are considerably larger than the DGLAP-based 
predictions particularly at low values of Q 2 . This demonstrates the need to better understand 
the role of higher twist (and power) corrections in parton analyses. Diffractive vector meson 
production is clearly a good process in which to further investigate these effects. 

In the parton-hadron duality approach we have a common mechanism for the description 
of all vector meson production processes, YP ~^ ^p, governed by the average hard scale (K 2 ), 
where 

K 2 = z(l-z)(Q 2 + M 2 ). (22) 

Therefore it is informative to plot all the observed cross sections, in a given W domain, as a 
function of Q 2 + My on the same plot, after allowing for the different photon-quark couplings e q 
(that is p : J/ip : T = 9 : 8 : 2) and the different energies W of the data. The result is compared 
with the p prediction in Fig. 3. The fact that the measured cross sections approximately lie on 
a common curve, demonstrates the universality inherent in the perturbative QCD description. 
Some departure from universality may arise from the different measured t slopes, from the 
flavour symmetry breaking of the qq — > V transitionf], and from comparing data with different 
average W values. 

Figs. 4, 5 show the QCD predictions for Or/or- The upper (lower) plot of Fig. 4 compare 
the ratio for p (J/i[>) electroproduction with the recent HERA data as a function of Q 2 at fixed 
energy W = 75 GeV (W = 90 GeV), whereas Fig. 5 shows the W dependence for fixed values 
of Q 2 for both p and J/ip production. Recall that in Ref. []TTJ it was pointed out that only 
0"^ is calculable in perturbative QCD; the calculation of 0r(p) using the leading twist p meson 
wave function is infrared divergent. We must therefore explain how the olI®t curves can be 
obtained? In the parton-hadron duality approach, with the J p = 1~ projection of the qq system, 
the integral over the quark kx is of logarithmic form Jl3| (like in a usual DIS amplitude). So 
the corresponding Feynman graphs have (at leading order) a pure ladder structure with strong 
kr ordering along the ladder. The factorization theorem is therefore valid for cry, as well as o^. 
After convolution with the gluon distribution, the logarithmic behaviour effectively enhances 
the transverse amplitude by a factor I/7, so ot ~ I/7 2 as in (|TJ). The decrease of 7 with 
increasing Q 2 masks the naive Q 2 /M 2 expectation for the Q 2 behaviour of (Jl/ot- 

6 In p4j the value of the unintegrated gluon is determined down to It — ko = 1 GeV. Below fco we use the 
linear approximation as was described in footnote 4 above, with xg(x, ko) = 1.57 (1 — x) 2,5 . 

7 In the non-relativistic approximation the qq — > V vertices are proportional to the meson wave functions 
evaluated at the origin, which may differ according to the mass of the meson. In our approach this is replaced 
by the different z intervals sampled by the relativistic p system as compared to the more non-relativistic J/t/j 
and T systems and by possible different choices of the mass interval AM covering the resonance peaks. Here 
the same value AM = 200 MeV was used for p and J/ijj production. For the T the interval M = 8.9 . . . 10.9 
GeV was chosen to predict T(15, 25, 3S) production in accordance with the experimental analysis and the 
resulting cross sections where divided by 1.7 to get the predictions for T(1S), see [ gjj . 
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It is interesting to compare the predictions for the W behaviour of <7 L /a T for J /if; production 
with those for p production, shown respectively by the dashed and continuous curves in Fig. 5. 
For p production we are in a relativistic qq situation where z covers an extensive part of the 
(0,1) interval allowing the I/7 behaviour to develop. The growth of o~l/o"t with W reflects the 
rise of 7 with 1/x. On the other hand J/ip production is nearer the non-relativistic limit where 
z = 1/2 and o~lIo~t = Q 2 /M 2 apply, and hence the ratio ctl/&t is almost independent of W. 

To gain physical insight we have discussed the results in terms of an effective gluon anoma- 
lous dimension 7 and the simplified formula (D). In the actual computations we use, of course, 
the explicit unintegrated (skewed) gluon distributions and perform the full integrations over kx 
and £t, or related variables. 

In summary we have shown that a perturbative QCD parton-hadron duality approach is able 
to describe all the main features of the <jl, and even the ctt, data for diffractive vector meson 
production j*P VVi provided that there is a sufficiently hard scale (that is provided Q 2 + M 2 
is greater than about 5 GeV 2 ). We emphasize that the approach is infrared convergent. There 
are non-negligible contributions at low gluon transverse momentum £t, but the perturbative 
gluon form matches well on to the linear l\ form for < £0 making the predictions rather 
insensitive to the choice of £ . The effects of using skewed gluons are fully included in the QCD 
calculations. The skewed distribution is completely determined by the conventional (diagonal) 
gluon distribution, and is found to enhance the p and J/ip cross sections by about 55% and 
70% respectively at the largest observed values of Q 2 . For p production the use of the skewed 
gluon distribution predicts a flatter Q 2 dependence, compatible with the recent data, see Fig. 2. 
We conclude that the data for diffractive vector meson production processes at HERA offer a 
particularly sensitive probe of the properties of the gluon distribution of the proton. The M, 
Q 2 , t dependences and the spin properties can be measured with increased precision and all 
constrain the behaviour of the gluon. 
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Figure 2: The QCD predictions for the Q 2 dependence of the cross sections for 7*p — > pp 
(upper plot, W = 75 GeV) and 7*p — > J/ijip (lower plot, W = 90 GeV) compared with the 
HERA data ||, |!|. The continuous curves are obtained using the MRST(99) gluon p2| . For 



the lower curve the default value 2K 2 is chosen for the scale of as in the K, factor, whereas 
for the upper curve the scale K 2 was used. The dash-dotted (dotted) curves show the results 
if the CTEQ(5M) |23| (KMS pi| ) gluon are used. The dashed curves show our results using 
the MRST(99) gluon and default parameters but without the effect of skewing. All predictions 
contain contributions from the real part of the amplitude as discussed in the text. The data 
point for J/ip photoproduction in the lower plot is interpolated between HI data for different 
values of W [@ and agrees well with the ZEUS result . 
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Figure 3: The data ||, ||, |j, |3], [7] for the 7*p — > Vp cross sections with V — p (circles), J/ip 
(triangles) and Y(IS') (squares: solid ZEUS, open HI, both slightly displaced from Q 2 = 
for readability) plotted versus Q 2 + My. The QCD predictions (with standard parameters 
as described in the text) are shown for comparison as continuous, dashed and dotted lines, 
respectively. The J/ip and T data (and errors) are corrected for (i) the different photon- 
quark couplings by multiplying the J/ip and T measurements by 9/8 and 9/2 respectively, (ii) 
the different W values according to the QCD predicted energy behaviour a(J/ip) ~ W lA (in 
agreement with the experimental measurements from ||) and cr(Y) ~ W 1 ' 3 . The upper dotted 
curve is obtained using a fixed slope parameter 6 = 4 GeV -2 , whereas the lower curve contains 
the slope as given in Eq. (p0|). 
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Figure 4: The upper plot shows the QCD predictions for the Q 2 dependence of <JlI<Jt for p 
electroproduction (at W = 75 GeV) compared with HERA data ||, 0], , partially at slightly 
different (average) values of W as indicated on the plot. The ZEUS measurement displayed by 
the open triangle is the one obtained by relaxing the s-channel helicity conservation condition, 
see |J. The different linestyles for the different gluons are chosen as in Fig. 2. Here the steeper 
continuous curve corresponds to the standard choice of 2K 2 as scale of as in the K, factor, the 
less steeper one to K 2 . Also displayed is the naive expectation ol/ctt = Q 2 /M 2 (steep dashed 
line). The lower plot shows o L jo T for J/ip production (at W = 90 GeV) compared to data 
from [|, §]. 
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Figure 5: The W behaviour of (y L ja T for fixed values of Q 2 for both p electroproduction 
(continuous curves) and J/ip electroproduction (dashed curves), obtained with our default 
parameters and the MRST(99) gluon. 
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